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Abstract 

Standard evolutionary theories of aging and mortality, implicitly based on assumptions of spatial averaging, 
hold that natural selection cannot favor shorter lifespan without direct compensating benefit to individual 
reproductive success. Here we show that both theory and phenomenology are consistent with programmed 
death. Spatial evolutionary models show that self-limited lifespan robustly results in long-term benefit to a 
lineage; longer-lived variants may have a reproductive advantage for many generations, but shorter lifespan 
ultimately confers long-term reproductive advantage through environmental feedback acting on much longer 
time scales. Numerous model variations produce the same qualitative result, demonstrating insensitivity 
to detailed assumptions; the key conditions under which self-limited lifespan is favored are spatial extent 
and locally exhaustible resources. Numerous empirical observations can parsimoniously be explained in 
terms of long-term selective advantage for intrinsic mortality. Classically anomalous empirical data on 
natural lifespans and intrinsic mortality, including observations of longer lifespan associated with increased 
predation, and evidence of programmed death in both unicellular and multicellular organisms, are consistent 
with specific model predictions. The generic nature of the spatial model conditions under which intrinsic 
mortality is favored suggests a firm theoretical basis for the idea that evolution can quite generally select for 
shorter lifespan directly. 


Author Summary 

Classic evolutionary theories hold that selection cannot directly favor reduced lifespan. Intuitively, selection 
should act against a gene contributing to the death of its owner. Mortality is instead explained as an indirect 
effect of traits with net benefit to the individual. However, these theories are implicitly based on averaging 
over populations. Spatial systems, in which heterogeneity can arise, have qualitatively different behavior. 
We show that in spatial population models, self-limited lifespan is favored even without direct compensating 
benefit, and that this phenomenon can help explain many observations in nature. In the models, longer 
lifespan is favored on short time scales, but environmental feedback based on heterogeneous resource use 
results in shorter-lived strains having a long-term advantage. This result can help to explain observations 
that have been surprising or difficult to explain with classic theories. These results add a key mechanism to 
our understanding of the evolution of lifespan. 

* j ust in. werfel@ wyss. harvard. edu 


1 



Introduction 


Lifespan extension is a topic of longstanding broad interest, with a long history of varied treatments advocated 
for those interested in living longer. Lifespans vary widely among species [IJ|^ and can be experimentally 
altered by genetic modification [5 -13 and selective breeding 14 15 . Standard evolutionary theory explains 


senescence as the result of two core effects, both based on the observation that the strength of selection 
decreases as age increases, since sources of extrinsic mortality reduce the older population size on which 
selection can act. The first effect, mutation accumulation 16 17 , describes traits with late-acting negative 


effects being only weakly selected against and hence not eliminated. The second, antagonistic pleiotropy 18 


and its disposable soma formulation 19 , observes that traits conveying both early-life benefits and later-life 
penalties can be actively selected for. These factors are generally considered to provide the complete basis for 
senescence and intrinsic mortality, while older ideas of evolution of direct lifespan control are considered 
untenable; there is a broadly accepted understanding that selection does not and cannot act directly in favor 
of shortening or otherwise limiting lifespan 20 -22 . This understanding accords with intuition: a gene that 


contributes to or hastens the death of its owner should be eliminated. 

However, a variety of observed phenomena have been and/or remain difficult to reconcile with accepted 
These include large single-gene effects on lifespan and other central regulators of aging 


frameworks 23 


semelparity ^e-time reproductffi^pparently triggering death, as in salmon ) [2pp , anomalous 


patterns of longevity JSL and senescence including apparently nonaging species m 27[|28| and even 


negative senescence 26,29 , aging in unicellular organisms 30 -33 , and significant plasticity of aging and 


lifespan 14 15 . At the same time, more recent work in spatial modeling of evolutionary systems has given 
new insight into altruistic behaviors of individuals within populations [34| , with relevance to phenomena 

social communication 


including reproductive restraint 35 36 


37 , and sexual reproduction 38 39 . Such 


models may thus help cast new light on empirical observations related to aging. 
In this paper we present and extend recent theoretical analysis 


40 showing an adaptive mechanism 


for explicit control of lifespan using spatial evolutionary models, and discuss the application of the results 
to empirical observations. In the models, heterogeneity of limiting resources and self-organizing popula¬ 
tion structures lead to an adaptive self-limitation of lifespan, without direct compensating benefit to the 
individual. The mechanism is different from the standard evolutionary genetic theories founded on muta¬ 
tion accumulation and antagonistic pleiotropy. In contrast to these classic mechanisms, which are based 
on selection directly on the reproductive success of individual organisms, this mechanism is based in spatial 
heterogeneity and it predicts the evolution of traits that may become advantageous only for descendants after 
many generations. To make this counterintuitive result clear, we show explicitly (Fig. that longer-lived 
mutants in prototypical simulations have a temporary advantage that may last for hundreds of generations, 
but eventually they are out-competed by the shorter-lived types. Thus, there is no compensating reproduc¬ 
tive benefit to the individual of the shortened lifespan, as there would have to be in a mean-field model. 
Our results show that the benefit that makes shorter life favorable need not vest through reproductive suc¬ 
cess of the individual, but rather through far descendants; moreover, the mechanism is generic in spatial 
environments. 

Why should spatial models be different from non-spatial ones? The key lies in the short-term advantage 
but long-term disadvantage incurred by overly exploitative variants 35 37 . “Selfish” variants that outcom- 


pete their neighbors, at the expense of depleting their local environment, have an advantage on short spatial 
and temporal scales. However, on longer scales, their descendants are left in impoverished environments 
and can be outcompeted by others in richer areas. A longer-range view of “time-dependent htness” 35 


recognizes these effects, and predicts the broad success of altruism and individual restraint in spatial mod¬ 
els. Non-spatial models, including age-structured population models, examine the evolution of mutations 
affecting lifespan in large, homogenous populations. They have not studied the effects of spatial resource 
heterogeneity, i.e., the distinction between persistent impoverished vs. rich environments, and thus produce 
qualitatively different results from those presented here. 

The numerical experiments we discuss here were designed to explore whether and how intrinsic mortality, 
a particularly extreme form of individual restraint, might be evolutionarily favored in spatial models. The 
results indeed showed that intrinsic mortality was consistently favored, for a wide range of model assumptions 
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Figure 1. Longer lifespan can give short-term advantage and long-term disadvantage in spatial 
models. Plots show relative reproductive success (number of descendants) of individual mutants with 
lifespan L, introduced into a steady-state population with equilibrium lifespan L = Lq (Fig. [^), as a 
function of time (left) and of generation number (right) in simulations of a consumer-resource model (see 
text). Parameter values are (A) g = 0.13, u = 0.2, Lq = 2.7, (B) g = 0.17, u = 0.1, Lq = 7.6, (C) g = 0.2, v = 
0.005, Lq = 170, with resource growth rate g, consumption rate u, reproduction cost c = 0, and no mutation 
{fj,p = fiq = 0) after the initial introduction. Error bars show the standard error of the mean among ten 
independent runs, each recording the mean for a set of 40,000 invasions. 
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Figure 2. Model snapshots showing different spatial distributions of consumers and resources. 

Resources are shown in yellow; consumers with intrinsic mortality (nonzero q) are shown in blue, those 
without it {q fixed at 0) in red; empty spaces are shown in black. The values of p and q for each panel are 
those which evolve in ascendance studies (Fig. 3B,C). Each panel shows 250 x 250 sites; parameter values 
for which populations without intrinsic mortality were not found to be stable are shown as empty lattices 
(details and full lattice sizes in the Models section). 


and parameters, with a very strong advantage over the lack of intrinsic mortality. Moreover, the results are 
consistent with a variety of empirical observations, including some which have posed challenges for traditional 
theories of aging and mortality. 


Results 

We considered a family of explicitly spatial models of an interacting pair of organism types, in which one 
population lives at the expense of the second, as in predator-prey, pathogen-host, or herbivore-plant systems. 
Since analytic treatments [40| average over or otherwise approximate the complex dynamic patterns that are 
responsible for the distinct behavior of spatial models, we relied primarily on direct simulation to characterize 
their behavior. The model (Fig. used a stochastic cellular automaton to consider a population of abstract 
organisms (“consumers”) with limited but self-renewing resources (e.g., prey or hosts). Sites in a two- 
dimensional lattice represented empty space, available resources, or a consumer together with resources (a 
consumer could not occupy a site in the absence of resources). At each time step, resources alone reproduced 
into neighboring empty sites, with probability g for each empty site; consumers reproduced into neighboring 
resource-only sites, with probability p for each resource site; consumers exhausted resources in their site, 
leaving empty space, with probability v per time step and c per reproduction; and consumers died due to 
intrinsic mortality, leaving resources, with probability q. A consumer’s mean intrinsic lifespan is L = 1/q. 
To explore the evolution of lifespan in the consumer population, the intrinsic death rate q as well as the 
reproduction rate p were heritable and independently subject to mutation: with probability pq {pp), a 
consumer offspring has a value of q (p) differing from that of the parent by ±eg (icp). 

Two types of simulations investigated the evolution of lifespan control in a population. In “ascendance” 
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Figure 3. Ascendance studies favor self-limited lifespan. (A) History of evolving consumer lifespan 
in one example simulation (g = v = 0.1,c = 0), showing population mean/maximum/minimum. (B,C) 
Steady-state average values of (B) consumer reproduction probability p and (C) intrinsic lifespan L = 1/q, 
for different values of parameters g and v and for populations with (solid) and without (dashed) intrinsic 
mortality. (D) Adding an additional source of extrinsic mortality of magnitude go decreases the evolved 
value of intrinsic mortality g, and increases evolved reproduction probability p. The unit of measure g' is 
the mean value of g that evolves for go = 0; g = n = 0.01, c = 0.4. All error bars show the standard error of 
the mean from ten independent trials. 


studies (Fig. [^, we tracked the evolution of the intrinsic death rate g and reproduction rate p over time 
in a randomly initialized consumer population, to investigate what values dominate in the long term. In 
“invasion” studies (Fig.|^, we considered the ability of a single consumer to take over a population, in order 
to ask: If a rare mutation could confer or remove the capacity for lifespan control, would that mutant have 
an advantage or a disadvantage in its later spread through the population? These studies introduced one 
consumer into a steady-state population, followed its lineage until fixation (extinction of either invaders or 
invaded), and examined the probability of successful invasion in 100,000 or more such trials. 

Additionally, in order to investigate the robustness of the results to model assumptions, we further per¬ 
formed ascendance studies with numerous variants of the model. These explored a variety of considerations 
relevant to the model’s applicability to real-world systems, to test whether the results were sensitive to these 
factors. Variants included: 

• explicit deterministic lifespan and 

• other age-dependent senescence patterns, rather than a fixed probability of death per time step; 

• local consumer mobility through rearrangement and/or migration; 

• increased dispersal range of consumers or resources; 

• resource depletion that is deterministic and gradual, rather than stochastic and binary; 
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Figure 4. A successful invasion of consumers without intrinsic mortality by those with the 
capacity for programmed death. (A) The fraction of invaders in the population (solid line) increases 
almost monotonically with time. The region dominated by consumers with the capacity for programmed 
death grows steadily: the dashed line shows the area of a circle (under periodic boundary conditions) whose 
radius increases at a constant rate (correlation r — 0.997). Resource growth g = 0.05, consumer consumption 
V = 0.2, consumer reproduction cost c = 0. (B) Snapshots at 150, 1250, and 2350 time steps (colors as in 
Fig.[^. 


• ability of consumers to adjust their rate of living in response to resource availability; 

• spontaneous resource generation; 

• consumers reproducing to replace neighbors in occupied sites; 

• resources reproducing while exploited by consumers; 

• sexual reproduction by consumers. 

The ascendance studies showed that self-limited lifespan was consistently favored, for all parameter values 
tested. Changing the values of the constant parameters g and v changed the population size and distribution 
of consumers and resources (Fig.[^ and specific equilibrium value of lifespan that evolved (Fig.|^); but the 
qualitative result, selection for lifespan limitation, was found in all cases. 

The capacity for intrinsic mortality changed the characteristic population structure, as qualitatively 
visible in Fig. Table shows quantitatively that the mean patch size (defined as the number of consumers 
in a connected set) was consistently larger for populations without intrinsic mortality than for those with it, 
for the same values of the ecological parameters {g,v}. 

Adding extrinsic mortality (corresponding to deaths through, e.g., predation or accidents) reduced the 
evolved level of intrinsic mortality q, such that the total mortality level favored by selection for a given set 
of parameters remained approximately constant (Fig. §3). 

The model variants showed that the key finding of selection for lifespan limitation was robust to changes in 
model details and assumptions. In almost every case, the outcome was the evolution of self-limited lifespan. 
The exceptions were consistent with this understanding and with previous work: When consumer dispersal 
was large compared to the size of the space, the system was effectively well-mixed and lifespan control did 
not evolve. When resource availability was so great as to support the maximum amount of consumption and 
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Table 1. Intrinsic mortality changes characteristic population structure. 



V = 0.005 

V = 0.01 

V = 0.05 

V = 0.1 

V = 0.2 

g = 0.005 

Immortal 

Mortal 

N/A 

8.7 ±0.4 

N/A 

7.2 ±0.4 

N/A 

1.01 ±0.01 

N/A 

2± 1 

N/A 

2.0 ±0.1 

g = 0.01 

Immortal 

Mortal 

25 ± 1 

8.2 ±0.3 

12.5 ±0.3 

7.9 ±0.3 

N/A 

1.01 ±0.01 

N/A 

1.02 ±0.01 

N/A 

2.0 ±0.5 

g = 0.05 

Immortal 

Mortal 

250 ± 50 
11.8 ±0.4 

79 ±8 
11.2 ±0.6 

12.3 ±0.5 

7.8 ±0.4 

8.2 ±0.7 
1.08 ±0.03 

6.3 ±0.4 
1.03 ±0.01 

5 = 0.1 

Immortal 

Mortal 

700 ± 300 

20 ± 1 

230 ± 50 

18 ± 1 

24 ± 1 

12.2 ±0.5 

11.6 ±0.9 

8.1 ±0.3 

7.4 ±0.7 
1.07 ±0.03 

5 = 0.2 

Immortal 

Mortal 

5000 ± 4000 
54 ±5 

700 ± 300 
46 ±6 

50 ±4 

28 ±2 

21 ±2 

17 ± 1 

10.5 ±0.6 

8.5 ±0.2 


Mean consumer patch size, averaged over all consumers, for different values of resource growth rate g and 
consumption rate v, and for populations with and without intrinsic mortality (“mortal” and “immortal”, 
respectively). 


reproduction the model allowed, the system evolved to these limits. In real-world systems, both dispersal 
and resources are normally limited. 

The invasion studies showed a very strong population advantage to the capacity for intrinsic mortality 
(Table . Consumers with this capacity (g = 0 for the invader, but potentially nonzero for its descendants 
through mutation) invading populations of consumers without intrinsic mortality (q fixed at 0 for all descen¬ 
dants) had a success rate typically two to three orders of magnitude greater than that of invaders without 
intrinsic mortality, across all values of the ecological parameters g and v tested. Conversely, consumers with 
no intrinsic mortality managed no successful invasions of populations with programmed death in a total of 
several million trials. 


Discussion 

The simulation results provide theoretical support for the idea that direct lifespan control, and programmed 
mortality and senescence as ways of achieving it, are consistent with natural selection. In contrast, it is 
widely reported that theory is incompatible with the evolution of explicit lifespan control 20 21 23 . This 
perspective has guided and constrained the interpretation of empirical findings. Our results suggest that 
classic mechanisms relevant to the evolution of lifespan are incomplete. The mutation accumulation and 
antagonistic pleiotropy theories provide established mechanisms contributing to gradual senescence; direct 
selection provides an additional mechanism, not a replacement—where multiple mechanisms are possible in 
biology, typically all may play a role. Our results suggest the broad applicability of an additional mech¬ 
anism missing from previous theories. This mechanism can provide additional interpretations and novel 
understanding of a variety of observations in nature. 


Empirical phenomena: providing needed explanations 

Direct selection for intrinsic mortality and senescence—not just selection for an individual benefit with 
senescence as a side effect—can be used to help understand empirical phenomena, particularly for cases that 


have posed problems for traditional theory 23 but are straightforward to explain with direct selection. Here 
we briefly discuss three such cases. 

(1) Organisms of some species reproduce only once and die after reproducing (semelparity), and may live 
for a variable period before the act of reproduction apparently triggers their death. Because programmed 
death is inconsistent with classic theories 23 , other interpretations for semelparity have been sought; the 
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Table 2. Consumers with the capacity for intrinsic mortality dominate those without it in 
invasion experiments. 


g = 0.05, w = 0.05 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.51 ±0.05) X 10-^ 

(1.0 ±0.7) X 10-'’ 

Immortal 

(2± 1) X 10-'> 

0 


g = 0.05, u = 0.1 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(1.91 ±0.04) X 10-^ 

(1.1 ±0.3) X 10-4 

Immortal 

(1.2 ±0.3) X 10"^ 

0 


g = 0.05, u = 0.2 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(1.23 ±0.04) X 10-^ 

(1.6 ±0.4) X 10"'* 

Immortal 

(2.8 ±0.5) X 10-^ 

0 


g = 0.1, u = 0.05 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.51 ±0.05) X 10-^ 

(2± 1) X 10-5 

Immortal 

(2± 1) X 10"^ 

0 


g = 0.1, V = 0.1 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.38 ±0.05) X 10-^ 

(5 ±2) X 10-5 

Immortal 

(1.4 ±0.4) X 10-^ 

0 


g = 0.1, V = 0.2 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(1.80 ±0.04) X 10-^ 

(1.1 ±0.3) X 10-4 

Immortal 

(1.1 ±0.3) X 10-^ 

0 


g = 0.2, V = 0.05 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.24 ±0.05) X 10-^ 

(2 ± 1) X 10-5 

Immortal 

(3 ± 2) X 10-5 

0 


g = 0.2, V = 0.1 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.44 ±0.05) X 10-^ 

(5 ± 2) X 10-5 

Immortal 

(1.1 ±0.3) X 10-^ 

0 


g = 0.2, V = 0.2 

Invaded 

Immortal 

Mortal 

Invader 

Mortal 

(2.12 ±0.5) X 10-2 

(1± 1) X 10-5 

Immortal 

(7 ±3) X 10-5 

0 


Probabilities of successful invasions for different values of resource growth rate g and consumption rate v 
(all with reproduction cost c=0), and for each combination of invaders and invaded having or lacking the 
capacity for intrinsic mortality (“mortal” and “immortal”, respectively). 



traditional explanation is that so much energy goes into reproduction that none is left over for essential 
physiological maintenance afterwards 23 . However, cases such as octopus where removal of the optic gland 


allows continued survival and resumption of copulation 41 are difficult to reconcile with this interpretation, 
and are simpler to explain with an explicit programmed death mechanism. 

(2) It is widely reported as a key prediction of traditional theories that higher extrinsic mortality rates 
should result in the evolution of shorter intrinsic lifespans 18,20,23,25 42 -44 ; others hold that extrinsic 


mortality should have no effect on evolved lifespan 45 . In contrast to both, a key empirical study observed 


that guppies that evolved subject to higher levels of predation exhibited longer lifespans and lower rates of 
aging 25 . The spatial model predicts this result because additional extrinsic mortality reduces the need for 
intrinsic mortality [^. Fig.[^ shows quantitatively that adding a source of extrinsic mortality unrelated to 
resource exhaustion (with probability qo), by replacing the intrinsic death probability q with {q + qo), leads 
to lower evolved intrinsic mortality. 

(3) It has been argued that in multicellular organisms—where the distinctions between gametes and 
other cells (germline versus soma), and between parent and offspring, are clear—senescence is inevitable, 
while unicellular organisms (lacking those distinctions) should be effectively immortal [l| [T^[^[^[47] . How¬ 
ever, certain multicellular organisms have been reported to have negligible or even negative senescence, i.e., 
displaying insignihcant or declining effects of aging according to population-statistical or physiological mea- 

[^[| [3l[^ [48l^ 


sures 


54 . Conversely, senescence has been reported in organisms with no clear germ/soma 


distinction [55] , and programmed cell death reported in unicellular organisms 30 -33 . The difficulty of 


reconciling such empirical observations with traditional theories has been pointed out 26 30 31 50 55 56 


These phenomena are straightforward to explain if mortality and senescence are determined by selection, 
with favored lifespan set according to ecological conditions regardless of whether organisms are single-celled 
or multicellular. The spatial model predicts that negligible senescence may be favored under appropriate 
ecological conditions (e.g., well-mixed populations, or those kept in check entirely by factors other than 
resource limitations); conversely, programmed death can be favored for unicellular consumers as easily as for 
multicellular ones. 

We note that this potential uncoupling of intrinsic mortality and multicellularity has implications for our 
understanding of evolutionary history. Programmed cell death, critical during development in multicellular 
organisms (morphogenesis), has traditionally been viewed as having multicellularity as a prerequisite—the 
sharp distinction between somatic and germ cells intuitively enables some of the former to die without 
compromising their genetic contribution to the next generation—and accordingly it has been argued to have 
arisen contemporaneously or later 47 57,58 . Our results suggest that programmed cell death likely arose 


first in unicellular populations [M ,59 -62 and later become co-opted for morphogenesis, perhaps helping to 
enable the evolutionary appearance of multicellularity. Similarly, the reproductive senescence observed in 
some unicellular organisms could have been an exaptation for the cessation of division in somatic cell 
lines. 


Empirical phenomena: providing alternate interpretations 

Finding lifespan control in a robust theoretical model opens the door to interpreting a wide range of ob¬ 
servations of aging-related phenomena as arising from natural selection. Such a direct-selection mechanism 
does not preclude classic mechanisms such as mutation accumulation and antagonistic pleiotropy, but rather 
augments them, providing additional explanations that may contribute to observed phenomena. Here we 
suggest examples of cases where future study may find that selection directly for lifespan control plays a 
significant role. Since the model conditions essential to the qualitative result of self-limited lifespan—spatial 
environment, exhaustible/renewable resource, and limited dispersal—are typical features of the real world, 
we would expect lifespan control to be a widely present phenomenon. 

(1) The great range of natural lifespans observed among organisms, both dissimilar ones and those 
otherwise similar |^|^, is traditionally explained as selection for longer life being insufficiently strong in 
the cases of shorter-lived organisms to extend their lifespans to match those of longer-lived ones. Instead 
or additionally, shorter lifespans could be actively selected for in some organisms, according to ecological 
conditions. 
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(2) Large single-gene or few-gene effects on lifespan were initially surprising from the viewpoint of 


traditional theories 23 


18 , where 


Some have been interpreted as examples of antagonistic pleiotropy 
increased lifespan carries a fitness cost in early life [^. A more direct interpretation is that they could 
represent effective means of control to regulate the tuning of lifespan for different circumstances. 

(3) The evolutionary loss of mouthpart function in many adult insects, such as mayflies, is often the 
primary factor limiting lifespan [^. The traditional explanation for this loss of function is evolutionary 
irrelevance after reproduction; instead, it could have come about through active selection as a means of 
limiting lifespan. Similar mechanisms may be present in other animals: for instance, elephants go through 
six sets of molars in their lifetime, and starvation after the final set wears out may be the leading cause of 
death for older elephants [^. 

(4) The model predicts that reproduction rate rises with increasing extrinsic mortality rate (Fig. §3). 
This is consistent with a field study of an island population of opossums, free from predators, which were 
observed to live longer and reproduce more slowly than mainland individuals 


63 


interpreted as evidence of increased senescence in the presence of increased predation 63,64 


These trends have been 
However, an 


alternate interpretation is that the shorter observed natural lifespan on the mainland is a direct result of the 
high predation {intrinsic lifespan is rarely observed there, as most opossums are killed by predators, with 
barely a quarter of females surviving to produce a second litter 63 ), and the reproduction rate is higher to 
compensate. 

(5) The model predicts the evolution of long lifespans for low consumption v relative to resource avail¬ 
ability (Fig. [^, 1^). This is consistent with observations of animals such as marine fishes, which have 


longer lifespans and slower metabolisms at greater depths 65 ; crocodilians, which go from weeks to years 
between meals and are long-lived species characterized as exhibiting negligible senescence 67 ; and cave 
animals, which, compared to surface dwellers, 
evolving significantly extended lifespan 


68 


appear to have adapted to conditions of food scarcity by 
Similarly, populations limited by factors other than resource 
availability are predicted to evolve longer lifespans (e.g., as with the guppies subject to higher predation and 
having longer intrinsic lifespans 25 ). 


(6) The model demonstrates that in many circumstances, a strain exhibiting increased reproduction 
along with increased death rate has an advantage over slower-reproducing, longer-lived individuals. This is 
consistent with studies of cancer where rapidly replicating cells also have increased rates of apoptosis 69 , 
suggesting an advantage over other possible cancerous types as well as normal cell populations. 

(7) It is sometimes raised as an objection to hypotheses of programmed aging that while many genes have 


been found which have large effects on lifespan, none has been found which eliminates aging altogether 23 
However, if self-limited lifespan has as powerful an adaptive advantage under as general conditions as our 
results suggest, it is to be expected that the genetic basis for it would be solidly and robustly established, with 
mechanisms for flexibly tuning lifespan to adjust to changing conditions but no simple way of eliminating 
the control altogether. 

(8) The model demonstrates an evolutionary advantage associated with reducing resource use through 
shortening lifespans. Other mechanisms for reducing resource use (particularly age-dependent) may have 
similar effects. For instance, metabolism typically declines with age [70| , which may allow greater longevity 
while leaving more resources for descendants. A related issue is that of post-reproductive care of indirect 
offspring in humans (the “grandmother hypothesis”). A strong relationship exists between post-reproductive 
lifespan in women and their number of grandchildren [71| . (Note that a form of time-dependent fitness, mea¬ 
suring reproductive output at later generations, is necessary to describe this effect: there is no relationship 
between a woman’s lifespan and her number of direct offspring [71| .) The benefits associated with grand¬ 
mother care may compete against the benefits of increased resource availability through early death that we 
have focused on. In organisms where such care is present, its effects oppose selection for shorter lifespan. 


with potential relevance to the notably long lifespans of humans and perhaps some other primates 72 -74 


Time-dependent fitness 

A useful characterization of the long-term behavior of evolutionary models is the average number of de¬ 


scendants of a given type after an amount of time T —the time-dependent fitness 35 36 . This measure 
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Figure 5. Mean number of descendants of a consumer with specified p and q, introduced into a steady-state 
population with g = 0.1,v = 0.1,p = 0.24, q = 0, pp = pq = 0. Without intrinsic mortality ((7 = 0), a higher 
reproduction probability gives a short-term advantage to a lineage, but on longer timescales the average 
number of surviving descendants plummets. Intrinsic mortality allows faster reproduction, giving a long¬ 
term advantage despite a disadvantage that may last hundreds or thousands of time steps {p = 0.6, q = 0.2), 
or even favoring a lineage on all time scales (p=1.0,g = 0.2). 


quantifies the long-term survival of the type including effects of changes caused by the type to its environ¬ 
ment. Simulations to obtain the time-dependent fitness are similar to invasion studies, but without mutation. 
Fig.[T] showed results within populations of mortal consumers. Examples of results within a population of 
immortal consumers are shown in Fig.j^ Immortal consumers with p at its evolutionarily stable value (here, 
0.24) have a persistent fitness approximately equal to 1. Higher reproduction probabilities (e.g., p = 0.6) 
give a significant short-term advantage, for hundreds of time steps, but subsequently the average number 
of surviving descendants falls to 0. Mortal consumers (with faster reproduction, but programmed death) 
can have the opposite pattern, with a short-term disadvantage, but a significant advantage in the long term 
(e.g., p = 0.6, g = 0.2). 

This result illustrates that immediate advantage or disadvantage does not necessarily determine long¬ 
term advantage—rather, processes operating over hundreds or thousands of generations can overcome more 
immediate considerations. Thus, kin selection alone 75 is insufficient to account for these phenomena of 


restraint (self-limited reproduction or lifespan): Kin selection is concerned with same-generation interactions, 
and is sometimes extended to immediate descendants; to characterize kin selection as relevant to descendants 
arbitrarily far in the future is to imply that the kin-selectionist argument is equivalent to evolution itself, in 
which case the concept provides no additional conceptual or mathematical benefit. 

In other cases, mortal consumers can have an advantage over immortals in both the short and long term 
(e.g., p = 1.0, g = 0.2). 

Looking at number of offspring as a function of generation number rather than time (Fig. right-hand 
plots) provides a connection to population growth rate in traditional theoretical treatments. If x is the 
generation number and y{x) is the number of invaders of that generation (y-value in the plot), then (/(I) 
is the average number of direct offspring of an initial invader; similarly, y{x l)/y{x) gives the average 
lifetime reproduction per consumer in later generations. Dividing this quantity by the consumer lifetime, 
l/{v + kc + q) (taking all causes of mortality into account), gives r, the average number of offspring per 
consumer per time step. The population growth rate is this birth rate minus the death rate, r— (u-|-fcc-|-g). 
However, as those plots show, an initial advantage in population growth rate does not necessarily characterize 
long-term reproductive success. The traditional single-generation measures of fitness and population growth 
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rate are not sufficient descriptors of long-term success in spatial systems. 


Classic perspectives on lifespan control 


Mainstream evolutionary theory considers it well-established that selection does not and cannot act in favor 
of decreased lifespan (unless more than compensated by a concomitant increase in early-life fitness, in the 


antagonistic pleiotropy framework 18 ). A few characteristic quotations illustrate this point: 


“The way evolution works makes it impossible for us to possess genes that are specifically designed to 
cause physiological decline with age or to control how long we live. Just as an automobile does not 
have a built-in plan for decline written in its blueprints, we do not possess genetic instructions that 


tell our bodies how to age or when to die.” 21 


“Evolutionary theory correctly asserts that aging is not an adaptive trait.” 22 


• “There is a widespread but erroneous tendency to regard aging as programmed.... There are powerful 
arguments why [such a program] should not exist.... The suggested benefits from aging are ones that 
serve the interests of the species or group. Whenever a benefit at the group level is assumed to supersede 
the contrary interests of the individual, any evolutionary hypothesis must confront the problem of 
‘cheating.’ Individuals in whom the aging program was inactivated by mutation would benefit from 
the sacrifice of others, while enjoying any fitness advantage that might accrue from immortality. Such 


mutations would therefore be expected to spread.” 23 


• “Any hypothetical ‘accelerated ageing gene’ would be disadvantageous to the individual. It is therefore 
difficult to see how genes for accelerated ageing could be maintained in stable equilibrium, as individuals 
in whom the genes were inactivated by mutation would enjoy a selection advantage.” 


20 


• “Aging, a decline in condition with increasing age apparent as a reduction in survival and reproductive 
output, is apparently disadvantageous for the individual. Why then, does it exist at all? Is there some 
hidden advantage to aging? Evolutionary biologists tackled these questions early on and concluded 
that aging does not have a function and exists only because natural selection is less powerful late in 
life... .These two basic tenets, that aging is due to a declining force of natural selection and is not 


adaptive, still form the conceptual foundation of the biology of aging.” 76 


• “There is a striking discrepancy between the diversity of theory on the evolution of senescence and 
its treatment in the literature. Empirical evaluations of the evolution of senescence focus almost 


exclusively on the classical theory, as do recent reviews.” 25 


These views can be traced back to Williams, who, in the article in which he proposed antagonistic 
pleiotropy, wrote: 

• “Natural selection should ordinarily proceed towards lengthening life, not shortening it. Such selection, 
at the individual level, could conceivably be countered by selection at the population level, if senescence 
somehow favored group-survival... .The efficacy of such selection depends upon a rather complicated 
series of assumptions ... A theory based on the simpler and more widely applicable principle of selection 
within a group would be preferable, unless the assumption of effective between-group selection proves 
to be necessary.” 


18 


In this way Williams argued that selection for shorter lifespan should not be considered if any alternative 
explanation exists. Such a view risks blinding evolutionary biology to an important explanatory process, 
particularly with increasing evidence that selection above the individual level is an important evolutionary 


force 35 -37 77 -80 
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Context and future work 


Previous models that predict selection for increased mortality account for it through intrinsic physiological 
tradeoffs: e.g., by putting energy into reproduction at the expense of maintenance, an individual may achieve 
earlier reproductive increases that more than compensate for lost later opportunities. In our model, the 
tradeoffs that ultimately result in greater net reproduction are imposed through an extrinsic mechanism of 
environmental feedback: dying earlier can leave resources that can most effectively be used by descendants, 
with the net reproductive benefits appearing generations later. 

The idea that shorter lifespans can be and are selected for directly goes back to at least 1870 [^. It was 
later rejected based on theoretical arguments that evolution of such a trait opposed to individual self-interest, 
like other altruistic behaviors, must require group selection, whose applicability should be accepted only as a 
last resort 18 . Analytic 40 and experimental 81 studies that take a mean-field or spatially mixed approach 


support the conclusions of the now-standard theories. However, such spatially averaged systems exhibit qual¬ 
itatively different behavior from systems which, like the real world, possess spatial extent. In particular, 
spatial systems routinely demonstrate altruistic behaviors 34-37 82 which are not evolutionarily stable 


in spatially averaged (mean-field) models or well-mixed laboratory populations 82 . Previous 

models, some spatial, have demonstrated that selection for self-limited lifespan is not a theoretical impos¬ 
sibility 85,86 , under assumptions such as continual introduction of highly advantageous mutations 87 


pre-existing senescence in the form of decreasing fecundity 88 89 or decreasing competitive fitness 
increasing age, frequent deadly epidemics among physically and genetically close relatives 
group selection among nearly-isolated subpopulations 


90 with 


91 , or explicit 


92 . However, a generally applicable mechanism for 


the active selection of lifespan control without such limiting assumptions has not been previously demon¬ 
strated. 

The robustness of our result that self-limited lifespan is favored, under many variations of model de¬ 
tails and parameter values, suggests that selection in favor of shorter lifespan and genetically programmed 
senescence may indeed be a quite general phenomenon. As such, it may have acted on the ancestors of 
human beings, with strong implications for human medicine. If aging is a functional adaptation, rather than 
a collection of inevitable breakdowns or genetic tradeoffs, then effective health and life extensions through 
dietary, pharmacological, or genetic interventions 64,93 -100 are likely to be possible, with potential for 


significant impact (e.g., altering two genes extends nematode lifespan fivefold [^). The effects of aging may 
one day be treated through manipulation of an underlying mechanism rather than as disparate symptoms. 
The fact that theoretical understanding of evolution can play a critical role in guiding health research should 
motivate a wider reevaluation of the evidence in relation to the theoretical frameworks. 


Models 

In the base model whose results are plotted in the figures, the dynamics comprise a sequence of synchronous 
updates of a two-dimensional spatial array of cells. At each time step, the following events simultaneously 
take place for each site in the array: 

• An empty site has probability 

Pe^r = 1 - (1 - 

of transitioning to a resource-only site, where g is the probability per time step of a resource-only 
site reproducing into a given neighboring empty site, and Mr is the number of resource-only sites 
among the empty site’s four nearest neighbors. This expression is equivalent to the statement that 
each neighboring resource-only site has an independent probability g of “seeding” the empty site. 

• A resource-only site has probability 


Nc 

PR^C = l-\{{l-Pi) 


I 
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of transitioning to a consumer site, where Me is the number of consumers among its four neighbors and 
the values pi are the corresponding consumer reproduction probabilities. This expression corresponds 
to each neighboring consumer site i having an independent probability pi of trying to reproduce into 
this site. If more than one consumer does so, one is chosen at random to be the parent. The offspring 
has p and q initialized equal to that of its parent; then with probability Pp {pq), that value is increased 
by tp (cq), or decreased by that amount with equal probability. 

• A consumer site has probability Pc^e = v + kc oi transitioning to an empty site (all resources in 
the site are consumed, resulting in consumer death), where k is the number of offspring the consumer 
produces in this time step, and probability Pc^r = <? of transitioning to a resource-only site (death 
due to intrinsic mortality). 

Once every 100 time steps, the maximum, minimum, and mean values of p and q in the consumer population 
are recorded. 

Simulations were performed on lattices of size 250 x 250 unless the consumer population was not stable 
(quickly going to extinction in such cases), particularly a problem for immortal populations. Accordingly, 
some results shown in Fig. were run on larger arrays, each on the smallest of {250 x 250, 500 x 500, 
750 X 750} for which a steady-state population of consumers could persist. Increasing the lattice size further 
does not change the steady-state values of p and q. For immortal consumers, we used: 

• 250 x 250: g = {0.05,0.1, 0.2}, u = {0.1,0.2} 

• 500 X 500: g = {0.05,0.1}, u = {0.005,0.01,0.05}; g = 0.2, u = {0.01,0.05} 

• 750 X 750: g = 0.01, u = {0.005,0.01}; g = 0.2, u = 0.005 

In other cases {g = 0.005, v = {0.005,0.01, 0.05,0.1,0.2}; g = 0.01, v = {0.05, 0.1,0.2}), immortal populations 
were found not to be stable even on arrays of 2000 x 2000 sites. Mortal populations were stable on 250 x 250 
arrays in all cases, with the exception oi g = 0.005, u = 0.2 which required a 750 x 750 array. 

In most of the simulations presented in this work, we focus on the limiting case of no cost of reproduc¬ 
tion (c = 0), due to the greater clarity of interpretation in that case: e.g., in immortal populations, the 
reproductive restraint that evolves cannot be attributed to an individual effort to conserve resources, since 
reproducing is free to the individual. Increasing c does not change the qualitative result that restraint is 
favored in the long term. We discuss the effect of nonzero values of c further in the section on mean-field 
analysis below. 

In ascendance studies, we initialized the lattice randomly with each site having a 55% chance of being 
empty, a 40% chance of having resources only, and a 5% chance of having a consumer, with p and q in the 
latter case randomly chosen from a uniform distribution between 0 and 1. 

Values of p and q are only meaningful within a finite range, [0,1] and [0 ,1 — v] respectively. (The 
latter expression arises since starvation and intrinsic death are taken to be mutually exclusive possibilities.) 
Accordingly, mutants with values outside these ranges are set to the boundary values. To ensure that 
this operation does not cause artifacts in the steady-state value of q or p, we performed simulations that 
progressively reduced the size of mutations Cp and Eg. Following an initial 100,000 steps to achieve steady 
state, an additional 100,000 steps were performed during which Cp and e, were halved every 10,000 steps, 
ending with a final 50,000 steps to obtain steady-state behavior (Fig.|^; the values reported in Fig. were 
based on this final period. The resulting minimum value of g, as well as the mean, was significantly above 0 
in all cases in Fig.[^ supporting the finding that finite lifespan is consistently favored. 

In invasion studies, we initialized the simulation with a steady-state configuration (following 200,000 time 
steps) of the “invaded” population. One consumer was then chosen at random and converted to an invader, 
with p unchanged, q set to 0, and pq set to 0.1275 or 0 according to whether the invader was mortal or 
immortal. If invaders succeeded in taking over the entire consumer population, the model was reset on the 
next time step to a steady-state configuration of the invaded variant. If the invader’s offspring became extinct, 
or the steady state was reset following a successful takeover, a new consumer was chosen at random and 
converted to an invader. The simulation continued for 100,000 introduction attempts, or 500,000 attempts 
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Figure 6. Traces showing mean (dark) and maximum/minimum (light) values of p (blue) and q (red) over 
the full course of an ascendance study, with initial transient period {t = 0-100000) and intermediate period 
during which Cp and €q are reduced {t = 100000-200000). Statistics reported in Fig. 3 are based on the final 
50000 time steps of each simulation. 


for immortals invading mortals (due to such an invasion’s low chance of success). In control studies in which 
both invader and invaded were mortal, q was left unchanged when an invader is introduced. Simulations 
for each {g, v) pair were conducted on arrays of the smallest size for which immortal consumer populations 
could persist, as specified above. Table gives probabilities of successful invasion for several values of g and 

V. 

A modified version of invasion studies was used to evaluate the long-term reproductive success of different 
types, as shown in Fig.(see also the discussion of “time-dependent fitness”). The simulation was initialized 
with a steady-state configuration for given values of parameters g, v, c; all consumers had the corresponding 
equilibrium values of p and g, and mutation was turned off (/ip = Pq = 0). One consumer was chosen at 
random and converted to an invader, and its offspring followed as in the regular invasion studies above, with 
q set to a new value for invaders and pq still 0. For the plots on the left side of Fig. evaluating number 
of descendants as a function of time, the initial invader was chosen from the entire consumer population. 
For the plots on the right-hand side of Fig. 1, evaluating number of descendants as a function of generation 
number, the initial invader was chosen from the set of newly-born consumers, so that all its offspring were 
identified and assigned the correct generation number. The latter plots allow comparison of the number 
of offspring after many generations for strains with different lifespans, and are constructed as follows: each 
consumer in the invading population has a generation counter, starting with 0 for the original single invader, 
and with each offspring’s counter being one higher than that of its parent. The plots show the number of 
offspring of a given generation number, no matter when in time they live. 

Variations on the base model 

In this section we briefly discuss how a number of qualitative changes to the base model affect its behavior. 
We consider: deterministic rather than stochastic limitations to lifespan (explicit programmed mortality, 
or rapid senescence); increasing mortality with age (gradual senescence); spatial mixing within a consumer 
population and migration of consumers within resource areas; increased dispersal for reproduction of both 
consumers and resources; deterministic resource use with continuous-valued resources, rather than binary 
resource state and stochastic consumption; individual consumers’ ability to adjust their “rate of living”; 
resources that regenerate spontaneously, rather than being spread only locally by existing resources; consumer 
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reproduction supplanting others already present; resources that continue to regenerate even when exploited 
by consumers; and sexual reproduction by consumers. 

All these model variants support the generality of the conclusion that finite lifespan is favored, when 
resources are limited and dispersion is local. For each variant, unless otherwise specified, we performed 10 
independent simulations on a 250 x 250 lattice with g = v = 0.05 and c = 0. Table shows collected results 
for various cases. 


Deterministic lifespan. In this variant, the genotype specifies intrinsic lifespan directly as a fixed length 
of time: L time steps after a consumer is born, it dies, regardless of remaining resources. Including some 
variability in the lifespan, by choosing it from some distribution with mean at L and a moderate variance, 
does not change the results. 


Increasing mortality with age. 

given by the Gompertz equation 
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In this variant, a consumer’s probability of death at each time step is 
m{t) = TOoe®*, where mg is a constant (chosen here to be 0.1), t is the 
number of time steps since the consumer’s birth, and g is a heritable value as in the base model. In this case 
q is not restricted by the simulation to be nonnegative; still, it evolves to a limited value. 

Note that for other conditions, this model variant can predict the evolution of different senescence pat¬ 
terns. For instance, if intrinsic mortality is set to be initially much higher than the equilibrium found in the 
base model (e.g., mg = 0.5), then negative senescence can evolve. 


Consumer migration. In this variant, we allow consumer mobility. For this purpose, we split each time 
step into three successive stages: first, reproduction of both resources and consumers; second, resource 
depletion and consumer death through both starvation and intrinsic mortality; and third, mobility, in which 
each consumer (asynchronously, in random order) is able to move. 

We tested three types of mobility: (a) a consumer trades places with a randomly chosen neighboring 
consumer (if any exist); (b) a consumer moves to a neighboring, unoccupied resource site if such exists; (c) a 
consumer exchanges places with a nearby randomly chosen non-empty site of either type. In all three cases 
self-limited lifespan is favored. Increasing the mobility distance to a limited extent (e.g., choosing the target 
site in version (c) from a 7 x 7 square centered at the consumer’s original position) does not change the 
qualitative result. Allowing consumers to move to a nearby resource site when they exhaust the resources 
in their own site, in order to avoid the starvation during that time step that would otherwise occur, likewise 
does not change the qualitative result. 


Increased consumer or resource dispersal. In this variant, consumers or resources are not limited 
to reproducing into only the four neighboring sites. For ease of implementation, a modified version of the 
model was used for these studies: Instead of the base model’s synchronous update where each site was 
simultaneously updated based on its value and those of its neighbors, an asynchronous update was used in 
which sites were successively chosen in random order for updates of the following form: empty sites remained 
empty; resource sites remained resources, and had a probability g of also reproducing by converting one empty 
site to a resource site; consumer sites had a probability v of becoming an empty site and q of becoming a 
resource site {v + q < 1), and separately a probability p of reproducing by converting one resource site 
to a consumer site. For such resource or consumer reproduction, the offspring site was chosen randomly 
from all sites of the appropriate type such that both the row and column indices differed by at most Rr 
or Rq (for resources and consumers, respectively) from those of the parent. For Rr = Rc = 1, this model 
produces the same qualitative results as the base model. (Note that these dispersal ranges are both effectively 
greater than those in the base model: reproduction is possible there only in the 4-neighborhood, here in the 
8-neighborhood.) 

Increasing the consumer dispersal range Rc can disrupt the local neighborhood relationships that make 
it possible for selection to favor restraint. Large enough Rc makes the system effectively well-mixed. In such 
cases consumption increases through unchecked selection for faster reproduction and longer lifespan, until 
the consumer population exhausts all available resources and goes extinct. The dispersal range above which 
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Table 3. Evolved trait values for selected model variants and conditions tested. 
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based on 10 independent experiments. See text for details of variants. 




this occurs is a function of the ecological parameters g, v, c as well as the size of the simulation space. This 
is because g, v, c, Rc affect the length scale of the characteristic population structure; if that length scale is 
too large compared to the space, the spatial nature of the model breaks down. Thus a given value of Rc can 
result in extinction on lattices of a given size but allow evolution of restraint and a sustainable consumer 
population on larger lattices. For instance, with g = v = 0.05 and c = 0, a consumer population with 
Rc = 1 consistently goes extinct on a 50 x 50 lattice (in 5800 ± 6300 time steps), but survives (i.e., persists 
for at least 250, 000 time steps, evolving intrinsic mortality) on a 100 x 100 lattice; Rc = 2 consistently gives 
extinction on a 200 x 200 lattice (in 60000 ± 37000 time steps) but survival on a 250 x 250 lattice; Rc = 3 
results in consistent extinction on a 250 x 250 lattice (in 5400±4900 time steps), survival in 3 of 10 trials on 
a 500 X 500 lattice (with extinction in 100000 ± 80000 time steps in the other trials), and consistent survival 
on a 750 x 750 lattice. 

Increasing the resource dispersal range Rr increases the extent to which consumers generally have 
resources available, which affects the extent to which intrinsic mortality is favored. For instance, with 
g = V = 0.05, increasing Rr to 5 results in enough resource availability that intrinsic mortality is not fa¬ 
vored; increasing the consumption rate u to 0.1 reduces resource availability and returns the model to the 
regime where intrinsic mortality is favored, and increasing v to 0.2 results in the evolution of still higher 
intrinsic mortality rates. 


Continuous-valued resources and deterministic consumption. In this variant, we treat resources 
as continuous-valued and consumption as occurring with a fixed rate. All sites are characterized by the 
quantity of resources they contain, with a value from 0 to 1. Consumers deplete resources by an amount 
V per time step, and c per reproduction. When resources in a site reach 0, the consumer there dies. If 
the consumer dies prematurely due to intrinsic mortality, residual resources remain for future exploitation. 
When an empty site is converted to one containing resources, the resource value there is set to 1. Partially 
depleted resources are left unreplenished. Replenishing partially depleted resources slowly or quickly over 
time (to a maximum of 1) does not change the qualitative result. 


Consumers can adjust rate of living in response to resource shortages. To consider the ability 
exhibited by many organisms to adjust rate of living in response to environmental conditions (e.g., dauer 
formation in C. elegans^ dietary restriction 23 101 ), we explored a variant based on the continuous-valued 
resource variant described above. In this variant, if a consumer is in a site with resource value below a 
threshold T, it adjusts its consumption u, reproduction p, and intrinsic mortality q all by a multiplicative 
constant k. Additionally, depleted resources are gradually renewed, by an amount D per time step, in both 
resource-only sites and those occupied by consumers. 


Spontaneously generated resources. In this variant, we consider spontaneous appearance of resources 
at any location in space. This model corresponds, e.g., to certain plant-herbivore systems, where a plant 
such as grass can be cropped down to its roots and regrow. Empty sites in this variant become resource sites 
with probability g at each time step independent of whether other resources are located nearby. Resources 
are generated much more readily than in the base model, and so a much lower value of g produces an overall 
level of resource production more comparable to that in the base model. For high values of g, the consumer 
population is not limited by resource availability and lifespan control does not evolve. 


Reproduction supplants existing consumers. In the base model, existing consumers prevent repro¬ 
duction by neighbors into the sites they themselves occupy, in effect perfectly defending their territories. 
The opposite extreme is for new challengers always to defeat older competitors. In this variant, we allow 
consumers to reproduce not only into sites with resources alone, but also those with consumers already 
present, replacing the previous consumer. 


Resources reproduce even when exploited by consumers. In this variant, we explored the conse¬ 
quences of allowing resources to reproduce even when consumers are present—i.e., resource sites reproduce 
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with probability g in the absence of consumers and g' in their presence (in the base model, (/' = 0). 

A nonzero g' increases the overall level of resource availability. Hence this variant requires reducing g 
and/or increasing v to obtain results quantitatively comparable to those of the base model. For example, for 
g = g' = 0.05, u = 0.23, the steady-state value of q is close to that of the base model tor g = v = 0.05 and 
for g = V = 0.23. For a given g and v, increasing g' leads to lower steady-state values of q (longer lifespan), 
just as increasing g for fixed v does in the base model (Fig. 3). 

It is possible in this variant, if the level of resource growth is high enough compared to v, for consumers to 
have little enough impact on the spatial distribution of resources that selection does not limit reproduction 
probability (for immortals) or lifespan (for mortals). For example, with g = v = 0.05, reproduction and 
lifespan limits occur for g' < g/3 and not for g' > g/2. 

Sexual reproduction by consumers. In this variant, when a consumer reproduces, a second consumer 
is randomly chosen to be the other parent of the offspring, with the values of p and q for the offspring being 
the average of those of the two parents plus mutation as in the base model. When the second parent is chosen 
from nearby (e.g., from a 7 x 7 region centered on the offspring), finite lifespan is favored, as is reproductive 
restraint: q evolves to a value significantly greater than 0, and p to a value less than 1. When the second 
parent is drawn from the entire consumer population (i.e., a form of global dispersal), no restraint evolves: 
q evolves to 0 and p to 1. (In this latter condition, the changes in q and p in the population occur over a 
much longer time scale than in the base model, and population extinction follows the reduction of Cp and eg 
described above. Delaying the start of that reduction allows q and p to reach those extreme values.) 
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